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Verfahren und Vorrichtung zur Digitalisieaing sowie zur Datenkomprimiemng 
analoger Signal? 



Die vorliegende Erfindung betriffl ein Verfahren sowie eine Vorrichtung zur 
Digitalisierung analoger Signale und ein Verfahren sowie eine Vorrichtung zur 
Datenlcomprimierung analoger Signale. 

Viele Bereiche der digitalen NachrichtenUbertragung haben das gemeinsame 
5 Ziel, ein zundchst analog vorliegendes Quellensignal (Messwerterfassung, 
Audiosignalverarbeitung in Studioqualit^, usw.) unter Ausnutzung der Vorteile 
der digitalen Obertragung zu einem Empf^nger zu transportieren, um es dort 
In Form eines analogen Ausgangssignals wiederzugeben. HierfUr muss das 
Analogsignal zum einen digitalisiert bzw. quantisiert warden, zum anderen 
10 muss eine entsprechende Datenkompression erfolgen, um das Signal bei 
m6glichst guter Qualitat Uber ein digitales Obertragungssystem mit begrenzter 
Rate zu ubertragen. Nach einer PCM-Encodierung warden Oblicherweise 
Verfahren zur Datenkompression angewendet. die Redundanz innerhalb des 
Quellensignals sowie Irrelevanz bzgl. spezieller Eigenschaften des 
15 Konsumenten des Signals ausnutzen. Die zahlreichen bereits existierenden 
Verfahren zur Digitalisierung analoger Quellensignale verfotgen teilweise 
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grundsdtziich verschiedene AnsStze, sie alle lessen sich Jedoch in zwei 
Hauptkategorien einteilen: 

I) Verfahren. bei denen die rekonstaiierte Wellenform die ursprOngliche 
Wellenform sehr nahe approximiert. d.h. es erfolgt Iceine Ausnutzung von 
5 In-elevanz. Wir venAfenden hier den Begriff ^wellenfonnerhaltend" anstelle des 
oft angewendeten Begriffs ^verlustlos" in Bezug auf die Wellenfonmcodierung. 
da die Digitalisierung eines analogen Signals bei einer begrenzten Datenrate 
prinzipiell nicht ^veriustlos** ist (d.h. unendliche Entropie einer kontinuierlichen 
Zufallsvariable). 

10 ii) Nicht wellenfomnerhaltende Verfahren. Sie sind beispielsweise im 
Audiobereich von entscheidender Bedeutung, da hier oftmats nur der 
subjektive Hfireindruck am Empfdngerausgang entscheidend ist (z.B. 
Ausnutzung psychoakustischer Maskierungseffekte). Hier spielen 
Wellenfomiveranderungen (Amplituden- und Phasenverzenxingen), die das 

15 Signal durch Quantisierung und Komprimierung erfdhrt, i.d.R. nur eine 
untergeordnete Rolle. OblichenA/eise fUhrt eine Signalverarbeitung unter 
Ausnutzung von Irrelevanz zu einer rekonstruierten Wellenform, die sich sehr 
stark von der ursprUnglichen Wellenfomi unterscheidet und zudem oft mit 
einer groHen Signalverzdgerung behaftet ist (z.B. aufgrund von 

20 Spektraltransformationen Oder dquivalenten blockbasierten Methoden). Die 
Obertragungsqualitat eines solchen Verfahrens ISsst sich nicht durch einen 
StOrabstand 10 logio(SNR), SNR: Signal to Noise Ratio, im klassischen Sinne 
messen; sie muss vielmehr durch aufwendige Hdrtests mit trainierten 
Personen enmittelt werden. FOr viele Anwendungsbereiche (z.B. in der 

25 Messtechnik, Aufeeichnung von Wellenformen fQr eine weitere 
Signalverarbeitung zu einem spateren Zeitpunkt, Echtzeit-SignalQbertragung 
unter Anwendung digitaler Modulationsverfahren, die keine nennenswerte 
Signalverzdgerung zulassen z.B. fur drahtlose digitale BUhnenmikrophone) 
sind derartlge Signalcodierungsverfahren jedoch gSnzlich unbrauchbar. 
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Es ist somit Aufgabe der vorliegenden Erfindung, ein Verfahren und eine 
Vorrichtung zum Digitalisleren sowie zum Komprlmieren analoger Signale mit 
verbesserter Qualitdt vorzusehen. 

Diese Aufgabe wird durch ein Verfahren zum Digitalisleren analoger Signale 
5 geman Anspruch 1. durch eine Vorrichtung zum Digitalisleren analoger 
Signale gerndK Anspruch 12, durch ein Verfahren zum Komprlmieren 
analoger Signale gemSQ Anspruch 8 und eine Vorrichtung zum Komprimieren 
analoger Signale gemSiH Anspruch 17 gelOst 

Somit wird ein Verfahren zur Digitalisierung analoger Quellensignale 
vorgesehen, bei dem eine Z>-dimenslonale spharisch logarithmische 
Quantisierung der analogen Quellensignale erfolgt. 

Gemaa einem Aspekt der voriiegenden Erfindung wird ebenfalls ein Verfahren 
sowie eine Vomchtung zum Komprimieren analoger Quellensignale 
vorgesehen. bei dem eine Digitalisierung analoger Quellensignale gemaR dem 
oben beschriebenen Verfahren erfolgt und wobei eine differenzielie 
Pulscodemodulation durchgefDhrt wird. 

Somit wird ein wellenfomierhaltendes Verfahren vorgesehen, welches die 
folgenden Anforderungen erfQIlt: a) geringe Datenrate durch Ausnutzung von 
gOnstigen Packungseigenschaflen mehrdlmensionaler Gitter 
(Vektorquantisierung) sowie Ausnutzung von AbhSngigkeiten in der Folge der 
Abtastwerte aus dem Quellensignal, b) extrem hoher Dynamikbereich, d.h. 
das SNR ist konstant uber einen sehr groHen Aussteuerbereich von z.B. 60 
dB und mehr bzgl. der Kurzzeitvarianz eines analogen Quellensignals, c) 
hoher segmentweiser objektiv messbarer StOrabstand fOr kurze Segmente 
von Abtastwerten, vorzugsweise im Sinne des VerhSitnisses von Varianz des 
Nutzsignals zum mittleren quadratischen Fehler, d) Unempfindlichkeit 
gegenUber speziellen Signalparametern wie z.B. segmentweise 
Wahrscheinlichkeitsdichtefunktion usw. und am wichtigsten e) EinfQhrung 
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einer extrem geringen Signalverzfigerung in der GrOdenordnung weniger (bis 
zu 10) Abtastperioden. 

Die Erfindung beruht auf dem Gedanlcen, ein Verfahren zur hochauflOsenden. 
wellenformerhaltenden Digitalisierung anatoger Signale vorzusehen, wobei die 
Ubiiche sl^alare logarithmtsclie Quantisierung auf mehrdimensionale 
Kugellcoordinaten Qbertragen wird und sich somit die hieraus resultierenden 
Vorteile wie z.B. ein Iconstanter StSrabstand Qber einen extrem holien 
Dynamikbereicli bei einem selir geringen Verlust in Bezug auf die Rate- 
Distortion-Theorie ergeben. Urn die im Queliensignai vorliandenen 
statistischen Abhangiglceiten fUr einen weiteren StOrabstandsgewinn 
ausnutzen zu icOnnen, wird die Differentieile Pulscodemodulation (DPCM) mit 
sphSrisch logaritlimischer Quantisierung kombiniert. Das resuitierende 
Verfahren erzielt eine wirksame Datenreduktion mit einem hohen Langzeit- 
sowie KurzzeitstOrabstand bei extrem geringer Signalverz5gerung. 

Weitere Ausgestaltungen der Erfindung sind Gegenstand der UnteransprQche. 

Die Erfindung wird nachfolgend unter Bezugnahme auf die beigefUgte 
Zeichnung ndher erl£iutert. 

Figurl zetgt eine Darsteliung der Proportionalitdt eines 

Kreisbogensegmentes zum Radius, 
Figur 2 zeigt ein Diagramm zur Veranschaulichung der Anzahl von 

Quantisierungszellen, die einem Radius zugeordnet sind, 

normalisiert auf die Zelienanzahl pro Absoiutbetrag eines 

Abtastwertes in skalarer Quantisierung, 
Figur 3 zeigt einen Graphen zur Veranscliaulicliung eines 

Stfirabstandverlustes bei verscliiedenen A - Werten, 
Figur 4 zeigt einen Graphen zur Veranschaulichung des 

Stdrabstandes gemSfi der vorliegenden Erfindung, 
Figur 5 zeigt ein Blockschaltbild eines DPCIVI-Encoder mit 

RQckwdrtspradiktion gemSK der vorliegenden Erfindung, 
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Figure zelgt einen Graphen zur Veranschaulichung des Signal- 

Rauschabstandes gemft& einem weiteren 

AusfOhaingsbeispiel der Erfindung, und 
Figur7 zeigt Graphen der Signalpegel und des Signal- 

Rauschabstandes gemdfi dem weiteren AusfOhrungsbeispiel 

der Erfindung. 

Das Ziel der logarithmischen Quantlsierung 1st eine hohe Dynamik des 
Quantisierers, also ein welter Bereich des mittleren Signalpegels. in welchem 
das SNR und somit der maximale relative Quantisierungsfehler 

V! 



konstant ist. Hierbei bezeichnet Aq, die Breite des i-ten 

Quantisierungsintervalles und r, den Jeweiligen Rekonstruktionswert. Diese 

Au^be fQhrt zu der logarithmischen Quantlsierung z.B. im Sinne des A- 
Laws gemaa K. Trdndle. R. Wel&. EinfUhrung in die Pulscodemodulatlon. 
15 Oldenboui^ Verlag. Munich. 1974. FOr den Bereich mittlerer Aussteuerung 
ergibt sich Ulr R » 1 folgender StSrabstand gerndf^ N. S. Jayant, P. Noll. 
Digital Coding of Wavefomns. Prentice-Hali. Englewood Cliffs, New Jersey. 
1984. 

101og,o(SNR) = R ' 6.02dB + lOlog.o . 0) 

(1 + In(i4}) 



20 wobei R die mittlere Rate der Quantisiemng (bit/Abtastwert) bezeichnet und 
A den Qblichen Parameter der logarithmischen Quantlsierung darstellt; d.h. 
der maximale relative Fehler ist fUr Abtastwerte mit einem Betrag > l/A 
(bezogen auf einen Signalberelch der Quantlsierung von -1 bis -i-l) konstant. 
Damit wird ab einem mittleren Signalpegel des Quellensignals von 



wo 2005/013492 



6 



PCT/EP2004/008217 



-5, := 20log,o(I/i4)dB (2) 

ein StOrabstand gemdQ (1) erreicht. Oer Wert B, charakterisiert damit den 
Dynamikbereich der logarithmischen Quantisierung. Innerhalb dieses 
Dynamikbereiches ist der Stdrabstand zudem gdnzlich unabhflngig von der 
Wahrscheinlichkeitsdichtefunktion (WDF) des Queliensignals, so dass das 
Verfahren universelt einsetzbar ist. Ailerdings stellt der Term 
101og,o(3/(l + \n{A)y) im Verglelch zu einer gleichmSiftlgen Quantisierung bei 
Vollaussteuerung durch ein gleichverteiltes Signal gewissermaHen einen 
"StOrabstands-Verlust" durch Kompandierung dar, der fUr einen grolten 
Dynamikbereicli zu bezatilen ist. 

Nachfolgend wird eine Quantisierung in Kugeikoordinaten In D Dimenslonen 
beschriebenen. 

Die Sphariscii logarltlimische Quantisierung gehfirt zur Famiile der 
Vektorquantisierungsverfaliren (von z.B. R. Gray, David L. Neuhoff. 
Quantization. IEEE Transactions on Information Theory, pp. 2325-2383, Oct. 
1998, l^anfred Herbert. Lattice-Quantisieaing von Sprach- und 
Sprachmodellsignalen. Ausgewdhlte Arbeiten Uber Nacliriclitensysteme, Nr. 
79. Editor H.-W. SchQBier, Eriangen 1991). Durch eine kugeidichte Packung 
in vielen Dimenslonen kdnnen betrSchtliche Stfirabstandsgewinne bei der 
Quantisiemng erzielt werden, auch wenn keine statistischen AbhSnglgkeiten 
innerhalb der Signale ausgenutzt werden kOnnen. Bei sphdrischer 
Quantisierung wird ein Vektor x :=(jf,,...,Xo)mlt D Abtastwerten in 
kartesischen Koordinaten durch Polarkoordinaten u•.= {^p^,...,q>^_^^r) 
ausgedrOckt. Die D - 1 Winkel g>, sowie der Radius r sind durch folgende 

Gleichungen gegeben 0': imaglnSre Einhelt, arg (•) : Argumentfunktion, welche 
den Winkel einer komplexen Zahl im Bogenmafi angibt): 
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= arg 



arg(ac, e[-«^, + ;r) 



(3) 
(4) 



/€{2 D-l}. 



/-I 



(5) 



Die Rekonstruktion der kartesischen Komponenten aus dem Vektor u erfolgt 
5 hieraus gema& 



X, = r bi,^'S\ti{ip^,y\ /€{/) 2} 

jc, = r -6, •cos(^,) = r-fto » 



(6) 
(7) 



mit den Radien b/der "Breitenkreise** der Einheitskugel (Radius 1): 
Vi = U (8) 



= i flcosCfl?;), i€{Z>-2,...,0}. (9) 



Damit ist ausgehend von ^^.^ eine aulwandsgOnstige rekursive Darstellung 

10 gefunden. Eine effiziente Implementieaing von Vor- und RUcktransfonnation 
ist mit Hilfe des sog. CORDIC-Algorithmus (welcher in Voider, Jack. The 
CORDIC Trigonometric Computing Technique. IRE Trans. Electronic 
Computing, Vol EC-8, pp. 330-334, Sept. 1959 beschrieben ist) mit geringer 
Komplexitat mOglich. 

1$ Nachfolgend wird die Spharisch logarithmische Quantisierung beschrieben: 
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Damit bei der Quantlslerung in Polarkoordinaten die Vorteile der 
logarttlimischen Quantisierung, d.h. die Unabhdngiglceit des SNR von der 
WDF des Queliensignals, weiterhin genutzt werden IcOnnen, wird nun der 
Betrag (Radius) iogarithmisch gemdH den Vorschriften des A -Law quantisiert. 
5 Die Winl<et werden indlviduell glelchmdllig quantisiert, Jewells mit Wahl der 
Quantisierungsinten/alle in Abhdngigl^eit der bereits quantisierten Winkel 

q>f, I € {/ + 1,/ + 2,...,D - 1} liOherer Ordnungen, wie sie auch fOr die 

Rekonstruktion verwendet werden. Somit wird eine sehr einfache 
Implementierung der Quantisierung und der Signalrekonstruktion erzielt, 
10 wobei durch die iterative Vorgehensweise gem^R (3) bis (9) eine schrittweise 
Verarbeitung der einzelnen Koordinaten wIe tm Fail skalarer Quantisierung 
bewalirt wird. 

Aufgrund der Proportionalitdt der Ldnge eines Kreisbogensegmentes zum 
Radius sind die Anforderungen an die logaritlimische Quantisierung, d.h. die 

15 Proportionalitdt der Grdl^e eines Quantisierungsintervails zum Signalwert, 
bereits fUr gleiclim^aige Quantisierung der Winkel erfOllt. Fig. 1 zeigt ein 
einfaches Beispiei fOr zwei Dimensionen D = 2. Der zweite Term in (1) muss 
somit nur fUr eine der insgesamt D Dimensionen bezahit werden, hierin liegt 
die wesentliche Ursache fOr die Gewinne mittels sphdrisch logarithmischer 

20 Quantisierung. 

Urn die Implementierung der Quantisierung in Kugelkoordinaten zu 
vereinfachen, gehen wir von ndherungsweise kubischen Quantisierungszellen 
aus, d.h. die Oberfiache einer O-dimensionalen Kugel mit Radius 1 wird durch 
ein Netz aus D - 1-dimensionalen (Hyper-)Wurfeln quantisiert. Wie 
25 nachfolgend gezeigt wird, ergibt sich durch diese suboptimale Quantisierung 
der Kugeloberfiache im Vergleich zu einem optimal dicht gepackten Gitter 
jedoch im Grenzwert O -» oolediglich ein Verlust von 
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also ein Ratenverlust von ca. 1/4 bit/Abtastwert, der zugunsten einfacher 
Implementierbarkeit hingenommen wind. 

Im fblgenden wird beschrieben. wie die Aufteilung der M'* (d.h. M :~ 2!^) pro 
Quantislerungsschritt zur VerfOgung stehenden Quantisierungsstufen auf die 
5 einzeinen Quantisierungsintervaile fQr den Radius und die Oberfldclie der O - 
dimensionaien Einheits(i^yper-)l(ugei in optimaler Weise vorgenommen wird. 

FQr den iogarithmischen Bereich 

der A -Law-Vorschrift gelte wie Qbiich die folgende Kompressoritennlinie 

Kr) = ro • c • In — + 1 1 mit c := - — \— (10) 
\ ''o J l + lnA 



10 



FQr deren Ableitung gilt: 



A bezeichnet den freien Parameter des A -Law und ro wird als noch naher zu 
beschreibender Normierungsfalctor fQr den Radius eingefOhrt. Mit Mo 
15 Quantisierungsintervalien fQr den Radius (O-te Komponente des Velctors u ) 
erglbt sicli fQr die Breite der Quantisierungszeiie in Richtung des Radius 
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Man beachte, dass Ar(r) im betrachteten Bereich unabhfingig von ro ist. 

Urn bei der Quantisierung in Polarkoordinaten ein fihnliches 
Obersteuerungsverhalten zu erhalten wie bei kartesischen Koordinaten, wird 
der Maximaiwert Tq fQr den Radius derart nonniert, dass die somit entstehende 
5 D-dimensionale Kugel das glelche Volumen besitzt wie ein O-dimensionaler 
WQrfel der Kantenlflnge 2 (entsprechend einem Quantisierungsbereich 

Xf € [- in jeder Dimension), d. h. V^^^^ = V^^^^^ : 




mit dem Volumen der 0-dimensionalen Einheitskugel gemSH J. H. 

10 Conway and N. J. A. Sloane. Sphere Packings, Lattices, and Groups. 
Springer-Verlag. 3rd edition. 

"''(Dm- 

Damit gilt: 

'b=-w ^0 >1 VDeIN (14) 

15 Um die Bestimmung des resuitierenden StSrabstandes zu verelnfachen. wird 
im Folgenden r = 1 gesetzt, d.h. es werden diejenigen Quantisierungszellen 
betrachtet, welche auf der Oberfiache einer Kugel mit Einheitsradius liegen. 
Dies ist ohne BeschrSnkung der Allgemeinheit mdglich. da das SNR aufgrund 
der logarithmischen Quantisierung im Bereicli 
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unabh^ngtg vom Radius ist. FOr die Breite der Quantisiemngszellen erhfllt 
man nun 



A:=-J— . (15) 
A/„-c 



5 Durch die bereits angesprochene gleichma&ige Winlcelquantisierung wird die 
Oberfldche der O-dimensionalen Eintceitskugel In Zeilen eingeteilt, welche 

jeweils die Form von D - 1-dimensionaien WQrfeIn besitzen. FQr die 
Oberfiache dieser D - 1-dimensionalen Kugel gilt gemad J. H. Conway and N, 
J. A. Sloane. Sphere Packings. Lattices, and Groups. Springer-Verlag, 3rd 
10 edition: 

S^Po^r'''' mit Po^D a^ . (16) 



Nachfolgend wird die Aufteilung der Quantisierungsstufen auf Radius und 
Winkel beschrieben. 

FOr eine faire Aufteilung der zur VerfOgung stehenden 

15 Quantisierungszellen auf Radius und Oberfldche der Einheitskugel gelten die 
folgenden Oberlegungen: im Hinblick auf die Approximation durcli kubische 
Quantisierungszellen. d.h. ungefShr gleiche Abmessungen der 
Quantisierungszellen in alien Dimensionen (Radius und 
Kreisbogensegmente), muss die Oberfldche der Einheitskugel = 1) in 

20 gleiche (D - f;-dimensionale WOrfel mit identischer Kantenldnge A unterteilt 
werden: 
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(17) 



Da A/ =s 2^ Intervalle pro Abtastwert zur VerfOgung stehen, muss zusdtzlich 



gefordert warden. 

5 Mit (15), (17) und (18) erhdit man fOr die Anzahl der fQr die Quantisierung des 
Radius zur VerfOgung stehenden Quantisierungsintervalle 



In Fig. 2 sind die Quotienten /(A//2) der Anzahl an Quantisierungszellen, 
die dem Radius gewidmet sind, normalisiert auf die Zellenanzahl pro 
10 Absolutbetrag eines Abtastwertes in skalarer Quantisierung (O = 1), 
gegenUber der Dimensionszahl fUr mehrere Beispiele des Parameters A der 
A -Law-Kompandierung aufgetragen. Die Werte von A sind derart gewShlt, 
dass die Anwendung der ^4 -Law-Kompandierung zu kleinsten 
Quantisierungsintervallen fUhrt, die um den Faktor 

15 2^,A/ieIN kleiner sind, ats dies bei gleichmsaiger Quantisierung bei 
identischer Gesamtintervallzahl der Fall ware: 



(18) 



1 1 



(19) 




bj'{r<.\IA) = 2-^--;^ ,d. h. 

M ^ 



(20) 
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2*- = 



A 



(21) 



Demzufolge ist die AuflOsung fOr kleine Werte von r um Ait bit erhOht. 

Fig. 2 zeigt. dass sphdrisch iogarithmische Quantisierung bis Qber den Faktor 
8 hinaus meiir Intervaile dem Radius zuordnet. ais eine skaiare Quantisierung 
s mit einer aquivalenten AuflOsung im Bereicli sehr kleiner Signatwerte. Daher 
wird ein Gewinn von 3 bit/Abtastwert Oder 18 dB fUr sphdrisch logarithmisclie 
Quantisierung in diesem Bereich erreicht 

Unter Venwendung von (15) eriiait man die Kantenldnge der Quantisienjngs- 
zellen auf der Oberfiache der Einheitskugel: 



Die Konstante c ist durch den gewdiiiten Dynamikbereich bestimmt. siehe (2) 
und (10). 

Die aktuelle gieiclima&ige Quantisierung der Winkelvariablen fOr gleiche 
Kreisbogensegmente auf der Einheitskugel erhdit nnan mit l-1ilfe der rekursiven 
15 Gleichungen (8) und (9). Wir haben M/ Intervalie fOr die Koorelinate zur 
VerfQgung: 
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(22) 




(23) 




(24) 
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for €{D-2,i)-3,...,2} 



2nb, 



(25) 



mit \x\ e IN rgrSfite ganze Zahl ^ ;c mit x e IR^ und entsprechend (8) 
und (9) fUr quantisierte Winkel . 



Es gilt zu beachten, dass die Anzahl Mi der Quantisierungsintervalle in der 
Dimension /eine Funktion der in den Dimensionen 1...mO- 1 ausgewflhiten 
Quantislerungszeile ist und daher ausgehend von M0.1 iterativ berechnet 
werden kann. Eine Berechnung im Voraus ist nicht mSglich. 



Die Zuordnung eines Indexes e {o,l,...,M'' -l} zu der aktuellen 
10 Quantisierungszelle sowie die Rekonstruktion kann durch Anwendung 
verketteter Lookup*Tabellen, wie sie aus dem Shell Mapping, z.B. gemaR 
Robert F. H. Fischer. Preceding and Signal Shaping for Digital Transmission, 
pp. 258-281, John Wiley & Sons. Inc., New York, 2002. ISBN 0471 22410 3 
bekannt sind, durchgefuhrt werden: 

15 Die Zuordnung eines Indexes zum Radius ist unabhdngig vom Zeilenindex auf 
der Einheitskugel, somit kann nunmehr das Problem der Indexzuordnung fUr 
die Oberfl^che der Einheitskugel betrachtet werden. 

Die Startnummer des Indexes ist 0, die maximale Endnummer des Indexes ist 
- 1, entsprechend der Anzahl an Quantisierungszellen auf der Oberfldche 

20 der Einheitskugel. Mo^i ist die Anzahl der Unterkugein der Dimensionalitat 
D - 1 und kann aus (23) berechnet werden. Zugunsten einer schnellen 
Implementierung wird eine Lookup-Tabelle erster Ordnung mit Mo-i ElntrSgen 
No,...,Nmd'\''\ venvendet. Hierbei bezeichnet A// den kieinsten Index alter Zellen, 
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die zu dem Men Quantisierungsintervall for den Winkel ^o-i SehOren: 

iVo = 0 (26) 

ie{0 (A/a.,-2)} 

5 Somit existieren Nm - N, Zetien fOr einen festgesetzten f-ten Wert von 0p.f . 
FQr jedes Nb wird eine Lookup-Tabelle zweiter Ordnung angelegt. welche die 
indizes enthdlt, die wiederum den ideinsten Index alier Zellen des /-ten 
Quantisierungslntervalls fOr den Winitel ^^-2 innerhalb der O - 1- 
dimensionalen Unterlcugel bezeichnen, auf die AZ/Sich bezieht. 

10 0,0=0 (27) 
/ e {0 iM^,_^-i)}, 



J e 



jo (^o-2(^i>-.='-f]-2]| 



Diese Vorgehensweise fuhrt zu O - 2 Ordnungen von verlcetteten Lool<up- 
15 Tabellen, urn die D - 1 Winlcel zu indlzieren (die Lool<up-Tabel!e der Ordnung 
D - 1 far den Winl<el ^, wQrde aufeinanderfolgende Zellenindizes enthalten 
und muss deshalb nicht tabeliarislert werden). 
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Nachfolgend wird das Quantisierungsgerdusch und der StOrabstand ndher 
betrachtet. 

Solange eine ausreichend grofte Anzahl an Quantisierungszellen in D 
Dimensionen vorliegt, darf die (Ibiiche Approximation eines gleichverteiiten 
5 Quantisiemngsfehlers innerhalb dieser kubischen Zellen angewendet warden, 
wobei jede Zelle durcli ihren Mittelpunlct reprdsentiert wird. FOr die Leistung 
des Quantisierungsgerdusches gilt ndlierungsweise 

^'D. (28) 

Die sphdrisch logarithmische Quantisierung erzwingt A(r)=:ATin alien 
10 Dimensionen fOr — ^ r ^ r^, dalier gilt 




siehe auch (22). Wie gewQnscht ist das SNR unabhflngig von der Varianz des 
Signals in diesem Bereich. 

15 Durcli Einsetzen von (1 3) und (1 6) erhdlt man: 



^P) = --47r ((^/2)!r c 



(31) 
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Betrachtet man (29). so kann F (D) ats Veriust bezUglich der Rate-Distortion- 
Schranke fUr i.l.d. gaua'sche Zufallsvariablen (6 dB-pro-Bit-Regel) Interpretiert 
warden. Dieser Zusammenhang 1st in Fig. 3 fOr verschiedene Werte von A 
dargestellt. 

5 Mit IHilfe der Stirlingschen Nfiherung 



berechnen, der den Veriust gegenOber der Rate-Distortion-Schranke 
10 darstellt. welcher durch die suboptimalen kubischen Quantisierungszellen 
(gegenUber (D - 1 )-dimensionalen Hyperkugeln fOr die Quantisierung der 
Oberfiache einer D*dimensionalen Einheitskugel und D c/^ verursacht 
wird. MIX anderen Worten: Durcli spharisch logarithmisclie Quantisierung 1st 
es m&glich, den in (1) bescliriebenen Veriust infolge der Kompandierung bis 
15 auf lediglich 1.53 dB wieder zu kompensieren. Somit ist es durch die Wahl 
des Parameters A und D theoretisch mdglich. einen beliebig groHen 
Dynamikbereicli (Bereicli konstanten StOrabstandes Qber dem mittieren 
Signalpegel) zu erzielen. ohne nennenswerte Verluste im maximal 
erzielbaren Stfirabstand in Kauf nehmen zu mUssen. Der asymptotische 
20 SNR-Verlust von 1,53 dB entspriclit einem Ratenverlust von 1/4 
bit/Abtastwert. der zugunsten einer extremen Reduzierung der Kompiexitat 
akzeptiert werden kann. 

Da im Sinne der Zielsetzung einer iogarithmischen Kompandierung bewusst 
von keiner Kenntnis einer Wahrscheinlichkeitsdichtefunktion fOr i.i.d. 
25 Quellensignalwerte ausgegangen wird, stellt aber wiederum die Rate- 




lasst sich der Grenzwert 



limW) = — = -l,53dB 



(32) 
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Distortion-Funktion fUr die gau&*8che Zufallsvariable die untere Schranke fOr 
die erreichbare Rate (bzw. Verzerrung) bei gegebener Verzerrung (bzw. Rate) 
dar (Berger's upper bound der Rate-Distortion-Function gemd& T. Berger. 
Lossy Source Coding. IEEE Transactions on Information Theory, pp. 2693- 
5 2723, Oct. 1998). Somit sind unter den hier gestellten Anforderungen und 
Einschrfinkungen Uber 1,53 dB bzw. 1/4 bit/Abtastwert hinausgehende 
Verbesserungen olinehin nicht mehr mOgiich. 

Fig. 4 zeigt als Beispiei den Abstand des StOrabstandes von der Rate- 
Distortion-Schranke R - 6dB abhdngig vom mittleren Signalpegel (10 logio 

10 (Varianz)) i.i.d. gau&verteiiter Zufallsvariablen be! A = 48270 fOr 
unterschiedliche Dimensionszahien der sphSirisch iogarithmischen 
Quantisierung. Diese Simuiationsergebnisse erfQIIen exakt die theoretische 
Analyse gemaR Fig. 3. Auf den ersten Blick scheinen die sehr gro&en Werte 
fUr A , wie sie in den Beispieien der Fig. 3 und 4 sowie in den nachfolgenden 

15 Beispieien verwendet werden, unreallstisch und niclit implementierbar. Es 
gilt hier jedoch zu beachten, dass gewfihnlich wesentlich mehr als M 
Inten/alle mittels (19) fOr die Quantisierung des Radius (Mo > M/2) venA^endet 
werden, siehe Fig. 2. Selbst fOr sehr niedrige Raten, z.B. R = 4 bit/Abtastwert 
existieren somit sehr feine Quantisierungslntervaile fQr den Radius. DarOber 

20 hinaus stellt (10) eine invertierbare Funktion fOr Jeden Wert von ^4 > 1 dar 
und ist somit in jedem Fail sehr gut fQr die Spezlfizierung eines bestimmten 
nichtgleichma&igen Quantisierers geeignet. Vorzugsweise soilte dieser 
Ansatz fOr echte Wellenformcodierung bei extrem geringer 
Signalverzdgerung nicht fUr Raten unter 3 bit/Abtastwert angewendet 

25 werden. 

DarOber hinaus wird in Fig. 4 deutlich, dass durch zwei Effekte mit 
wachsender Dimensionszahl bei konstantem Parameter A der 
Dynamikbereich stark ansteigt. Zum einen nimmt der Normierungsradius Tq zu, 
zum anderen wirkt die Begrenzung der Iogarithmischen Kompression nur auf 
30 eine Dimension, auf den Radius, wShrend fur O - 1 Dimensionen die 
Proportionaiitat der Ausdehnung der Quantisierungszelle zum Signalwert 
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(sBogentflnge zu Radius bei fixer Winkeldifferenz) auch fOr klelnste 
Signalwerte erhalten bleibt Somit erweitert sich der Dynamikbereich in D 
Dimensionen nflherungsweise auf 

« B,+20Iog,o(ro) + IOIog,o(Z)). (33) 

5 Zudem steigt aufgaind des MittelungsefTektes innerhalb von Jeweiis O Werten 
die Obersteueaingsfestigkeit, wodurch eine weitere VergrO&erung des 
Dynanr)lkberelches erzielt wird. Im Grenzwert (O ^ oo) ergibt sich ein unendlich 
groKer Dynamikbereich fOr Jeden beliebigen Wert von A . 

Die strukturbedingte VerzOgerung von sphflrisch logarithnnischer 
10 Quantisierung betrSgt, wie fUr Jede Vektorquantisierungsmethode, exakt D 
Abtastwerte. Wie aus Fig. 3 und 4 ersichtlich 1st, wird der grO&te Teil des 
mOglichen Gewinns bereits bei sehr kleinen Werten von D (bis 5) erzielt. 

Nachfolgend wird die Kombination von einer sphdrischen Quantisierung und 
DPCM beschrieben. 

15 Korrelationen zwischen den Abtastwerten q[k] werden mittels differentieller 
PCIVI (DPCM) effizient ausgenutzt, siehe Fig. 5. Ein Pradiktionsfehlersignal 
x[k] wird durch die Subtraktion von prddizierten Abtastwerten, die mittels eines 
linearen PrSdiktorfilters Hp(z) - aus rekonstruierten Abtastwerten q[k\ 
gewonnen werden, erzeugt. FQr einen idealen PrSdiktor zeigt die 

20 Prddiktionsfehlersequenz xfk] ein weiRes Leistungsdichtespektrum (PSD) 
sowie minimale Varianz. Oblichenveise ist dieses Prddiktorfilter adaptiv, urn 
sich an eine nichtstationSire Quelle anzupassen. Vorliegend werden nur 
Beispiele fQr ein festes Prfldiktorfilter vorgestellt, mit dem Zlel, die Interaktion 
mit spharisch logarithmischer Quantisierung mfiglichst einfach zu 

25 beschreiben. Beispielsweise bei Audiosignaien bleten sich bereits sehr kurze, 
feste Pradiktionsfehlerfilter (entworfen bzgl. eines Kompromiss-Kriterlums) 
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Oblicherwelse Gewinne von mehr als 18 dB bei einer Abtastfrequenz von 44,1 
kHz. Selbst segmentweise Gewinne, gemittelt Uber 6000 Abtastwerte 
(0,136s), von weniger als 15 dB lessen sich nur sehr selten beobachten, siehe 
hierzu Fig. 7 im Vergleich mit Fig. 3. Dartlber hinaus ist der Gewinn durch 
5 adaptive PrSdiktion begrenzt. falls die Signalverzfigerung auf einige wenige 
Abtastwerte beschrflnkt ist und merkliche Effekte, welche durch 
Aktualisiemngen der Koefflzienten hen^orgerufen werden, vermieden werden 
mQssen. Selbstverstdndlich kann das folgende Material direkt auf adaptive 
Prddiktion verallgemeinert werden. 

10 Es gilt zu beachten, dass ein SNR-Gewinn, ausgedrOokt durch eine 
Verringerung des mittleren quadratischen Fehlers zwischen originalen und 
rekonstruierten Abtastwerten, bei Anwendung der sogenannten 
ROckwartsprddiktion, wie in Fig. 5 dargestellt, mdglich ist. 

Bei logarithmischer Quantisierung ist das SNR des Prddiktionsfehlersignals 
15 x[k] unabhangig von seiner Varianz oder anders ausgedrOckt, die Leistung 
des Quantisierungsgerausches ist proportional zur Signalleistung. Somit kann 
der Pradiktionsgewinn, d.h. der Quotient der Varianzen von q[k] und x[k], 
direkt in einen SNR-Gewinn QberfQhrt werden. Damit ist logarithmische 
Quantisierung eine gute Wahl fOr DPCM. Hinzu kommt, dass durch DPCM - 
20 verglichen mit PCM - keine weitere systembedingte Signalverzdgerung erfolgt, 
da ein optimales Pradiktionsfehierfilter fOr maximalen Pradiktionsgewinn 
prinzipiell kausal und strikt minimalphastg ist (vgl. L. Pakula and S. Kay. 
Simple Proofs of the Minimum Phase Property of the Prediction Error Filter, 
IEEE Transactions on ASSP, vol. 31, 1983) und daher ohne strukturbedingte 
25 Verzdgerung invertiert werden kann. 

Nachfolgend wird das Gradientenabstiegsverfahren naher beschrieben. 

FOr die Anwendung von spharisch logarithmischer Quantisierung auf DPCM 
mit Ruckwartspradiktion gilt die gleiche Problematik wie fOr jede andere 
Vektorquantisierungsmethode: fOr die Berechnung des aktuelien 
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Prddiktionsfehlersignalwertes x[kj mOssen alio vorausgehenden 
rekonstruierten Abtastwerte q[k-i], i = l»2,... vorliegen, fOr einen hohen 
Prddiktionsgewinn sind insbesondere die unmittelbar vorausgehenden Werte 
(/ = 1^ / s 2) unverzichtbar. Diese Forderung ist leider nicht vereinbar mit der 
5 Quantisierung von BlOcken der Ldnge O Abtastwerte. 

Als LOsung fUr dieses Problem kann jede in der LIteratur en/vdhnte Methode 
zur Kombination von Vektorquantisterung mit DPCM herangezogen werden. 
Hier wird eine Methode» in Aniehnung an das Prinzip "analysis by synthesis", 
venvendet, welches belspielsweise von CELP-Wellenformcodierungs* 
10 verfahren geman N. S. Jayant» P. Noll. Digital Coding of Waveforms. Prentice- 
Hall, Englewood Cliffs, New Jersey, 1984 bekannt ist und kombinieren diesen 
Ansatz mit einem diskreten Gradientenabstiegsverfahren. 

Um den Widerspruch zwischen DPCM und Vektorquantisierung zu 
bewaitigen, wird die quadratische Euklidische Distanz zwischen Vektoren von 
15 Abtastwerten 

q[/] = iq[l . DUV • D + 1] q[l • D + £> - 1]) (34) 

und einem entsprechenden Rekonstruktionsvektor minimiert, l = \k/Dj 

An dieser Stelle gilt es zu beachten, dass neben der sphdrisch 
logarithmischen Quantisierung die Berechnung des entsprechenden 
20 Prddiktionsfehlersignals sowie die Inversion des Pradiktionsfehlerfilters in die 
Berechnung eines Paares 9,9 mit einbezogen werden muss, wobei auf 
vorausgehende Rekonstruktionsvektoren ^[/-m],m = 1,2,... zugegriffen wird, 
diese selbst aber unverSndert bleiben. Das Ziel des Algorithmus ist es, 
diejenige Quantisierungszelle fUr ;c[/] zu finden, fUr welche die Metrik 



(35) 
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minimiert wird. Urn einen geeigneten Startwert fOr den Algorithmus zu finden. 
schlagen wir vor. mit einer Vorwdrtsprddiktlon fUr die aktuellen D Abtastwerte 
zu beginnen, dies entspricht einer Deaktivierung der Kette aus ADC und DAC 
In Fig. 5 bzw. der direkten Einspelsung von 
5 q[k],k = l'DJ D'¥l,...J*D'^D-\ In das PrSdiktorfilter (in Fig. 5) anstelle 
von q[k]. Nun wird der resultlerende Vektor x sphdrisch logarithmisch 
quantislert und dabel ein Startvektor y erzeugt, aus dem wiederum durch 
gewOhnliche Inversion des Prddiktionsfehlerfilters (DPCM-Empfdngerstruktur) 
ein Vektor q folgt. Auf diese Weise kann die Metrikberechnung fOr einen 

10 gegebenen quantlsierten Vektor >;,erfolgen. Ausgehend von einem aktuellen 

Rekonstruktionsvektor werden nun die Rekonstruktionsvektoren aller 2D 

ndchsten Nachbarn yj^^^ in D Dimensionen ermittelt sowie die zugehdrigen 

Metriken gemdfl (35) berechnet und miteinander verglichen. Der Vektor mit 
der geringsten Metrik wird fur die nSchste Iteration verwendet, d.h. arg miny 

15 d^iq^qj^f^) erzeugt eine Aktualisierung von /, wobei q^ dem quantlsierten 

Vektor entspricht und durch Inversion des Pradiktionsfehlerfilters berechnet 
wird. Existiert kein solcher Vektor, so bricht der Algorithmus ab und gibt den 
Index des zu Ubertragenden Vektors y^ aus. Hier ist anzumerken, dass die 

Nachbarzellen ubiicherweise sehr verschiedene Indizes haben und es nicht 
20 trivial ist, diese Zellen zu identifizieren. Betrachten wir beispielsweise D = 3, so 
werden die Winkel des Azimuths (^, ) bei grofien Elevationen {ipjnahe 
±7tl2) grOber quantisiert, als bei kleinen Elevationen ((p^ nahe 0), weii 
letztere einen grfil&eren Breitenkreis auf der Kugeloberfldche besitzen. FOr 
eine schnelle Implementierung kOnnten die Indizes der Nachbarzellen in 
25 einem ROU fOr alle Zellen gespeichert werden, solange nicht zu grofi ist. 

Da lineare PrSdiktion mit in den Optimierungsprozess einbezogen ist, ist das 
resultlerende SNR manchmal grOf&er als es durch die Addition der Gewinne 
(in dB) aus sphdrisch logarithmischer Quantisierung und Prddiktion (DPCM) 
vorhergesagt werden kann. Besonders fUr den Fall kleiner Raten (z.B. R < 5 
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bit/Abtastwert) sind bemerkenswerte SNR-Gewinne beobachtbar, was 
nachfolgend und in Fig. 6 beschrieben wird. Ein analytlsches Ergebnis fOr das 
durch diesen Algorithmus erzielbare SNR ist den Autoren derzeit nicht 
belcannt. 

5 Simulationen zeigten, dass die im Mittei erforderliche Zahl an iterationen bei 
ca. 0.25 (pro D Abtastwerte) liegt, und bereits bei einer Besclirankung auf 
max. 3 Iterationen kelne signifikanten Verluste im Vergteich zu einem 
unbesctirdnkten Suchraum zu verzetchnen sind, so dass die Suche nach der 
optimalen Quantisierungszelle fUr niedrige Dimensionszahlen durchaus in 

10 Eclitzeit implementierbar ist. DarOber hinaus liefem bereits kleine Werte von D 
groUe Gewinne, sielie Fig. 3 und 6. 

Empfangsseitig ergibt sich keinerlei Anderung gegenOber Qbiichen DPCI\^- 
Verfaliren. An dieser Stelle soli darauf hingewiesen werden, dass die 
GesamtverzOgerung des Obertragungssystems ledlglich D Abtastwerte 
15 betragt und es sich somit fOr Obertragungen. welclie den Anforderungen einer 
extrem kurzen VerzOgerung entsprechen mOssen, hervonragend eignet. 

Naclfifolgend werden Varianten des Verfahrens naher beschrieben. 

Zur Beschleunigung der Suche nach der gQnstigsten Quantisierungszelle, die 
zu einer minimalen Verzenrung z.B. gemSli (35) fuhrt. sind prinzipiell alle 

20 Algorithmen der sog. Lattice*Decodlerung bzw. zur Auffindung eines 
Maximum-Likelihood-Codeworts bei der Kanalcodierung anwendbar, vgl. z.B. 
Erik Agrell. Thomas Eriksson, Alexander Vardy, Kenneth Zeger. Closest Point 
Search in Lattices. IEEE Transactions on Infomiation Theory, pp. 2201-2214, 
Aug. 2002 und die dort angegebenen Referenzen. Dabei ist vorzugsweise 

25 eine Transformation dieser Verfahren in Kugelkoordinaten vorzunehmen. 

Eine Variante ohne iterative Bestimmung der Quantisierungszelle ist dadurch 
gegeben, dass anstelle einer BerUcksichtigung statistischer Abhdngigkeiten 
innerhalb der aktuell zur Quantisierung anstehenden D Signalwerte durch ein 
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lineares Prddiktionsfilter gemflU Fig. 5 (also DPCM)» eine nichtgleichma&ige 
Quantisierung der Winkeikoordinaten vorgenommen wird. Dabei sind die 
linearen statlstischen Bindungen der aktuell zur Quantisieaing anstehenden D 
Signatwerte zu vergangenen Signalwerten mittels klassischer DPCM mit 
5 RQckwdrtspradiktion zu nutzen. Es entsteht also eine RUckwartspradiktion, bei 
der nach O Signalwerten eine Aktualisierung des Pradiktionsfilters urn O 
Schritte erfolgt. 

Es wird vorgeschlagen, die logarithmlsche Quantisierung des Betrags aufrecht 
zu erhalten, so dass anschlle&end die Signalpunkte so normiert werden 

10 kfinnen, dass sie auf einer Kugel mit Radius 1 liegen. Dies eriaubt die 
analytische Berechnung der Wahrschelnlichkeitsdichtefunktion der 
Signalpunkte auf der Kugeloberfiache, z.B. be! Annahme eines gauK'schen 
Signalprozesses aus der Autokorrelationsfunktion des Quellensignals, bzw. 
eine direkte experimentelle Ermittlung von relativen Haufigketten der 

15 Signalpunkte auf der Kugeloberfiache. 

FQr M » 1 lasst sich hieraus z.B. mit Hilfe des Optimierungsansatzes nach [9] 
eine nichtgleichmaRlge kubische Quantisierung der Oberfiache einer Kugel 
mit Radius 1 bestimmen. Dieser Ansatz besagt, dass im Mittel der Beitrag 
jeder Quantisierungszelle zum Quantisierungsgerausch gleich sein soil. 
20 Bezeichnet A; die Kantenldnge einer (O - 1)-dimensionalen kubischen 
Quantisierungszelle sowie z/ den dazugehOrigen Mittelpunkt, also den 
empfangsseitigen Rekonstruktionsvektor mit dem Betrag 1, sowie /(v);v := 

(^i,...,^0.,)die Wahrscheinlichkeitsdichtefunktion bzw. relative Haufigkeit von 
Signalwerten auf der Kugeloberfiache, so sol! also naherungsweise erfUllt 
25 sein: 
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(36) 



mil der Nebenbedingung 2^/"'* = 



(37) 



Hieraus fotgt unmittelbar eine (O - 1)*dimensionale Kompressor-Funktion 
^(v), durch welche die ungleichmaHige Quantisierung der KugeloberflSche 
5 eindeutig bestimmt wird. Die ungleichmdHige Quantisieaing kann dann 
beispielsweise wie im eindimensionalen Fall durch eine nichtlineare 
Verzerrung des Vektors v in einen Vektor z:-k{y), anschlie&ende 

gleichmaiSige Quantisierung zum Vektor gemail Abschnitt lll-C und 
Anwendung der Umkehrfunktion v := k'\z^) nachfolgend zur Gewinnung des 
10 Rekonstruktlonsvektors vollzogen werden. 

Vorzugsweise ist die resultierende (D - 1 ) -dimensionale Kompressorfunktion 
durch eine analytisch beschrelbbare Funktion anzundhern, Oder, analog zu 
der bei der eindimensionalen Kompresslon Qblichen Darstellung mittels 
GeradenstQcke (z.B. 13-Segment-Kennllnie, siehe [5]), durch (D - 2)- 
15 dimenslonale Teilebenen (mil konstanten partiellen Ableitungen) zu 
approximieren. Bei einer Qber dem Raum 



orthogonalen Rasterung dieser Approximation Qber den O - 1 Winkein kann 
eine einfache Auswertbarkeit der Kompressorfunktion erreicht werden. 

20 Zur Simulation von sphdrisch logarithmischer Quantisierung in Kombination 
mit DPCM wurde die Quantisierung der Ouverture sowie der Arie "Der 
Vogelfdnger bin ich ja" der Oper ^uberfldte" von Wolfgang Amadeus Mozart 
(vgl. Philips Classics Pruductions 1994 (DDD): Mozart: "Der Vogelfanger bin 
ich ja" (Die ZauberflOte). Polygram Records #442569-2, Track 3) untersucht. 
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Fig. 6 zeigt Simulationsergebnisse fOr f? = 3, /? = 4 und R - 7 bit/Abtastwert 
bei A= 102726. A = 48270, A = 4858 sowie verschiedenen 
Dimensionszahlen O. 

Das Signal wurde zuerst codiert. anschlieHend decodiert und das SNR durch 
5 Vergleich mit dem originalen CD-Signal berechnet. Der StOrabstand ist hier 
Jewells Qber das gesamte MusikstOck gemittelt. Bei den niedrigen Raten steigt 
das SNR mit D im Vergleich zu Fig. 3 noch steiter an, Ursache hierfOr ist eine 
verbesserte Mittelwertbildung. Beim Vergleich der Ergebnisse sind die 
verschiedenen Werte das Parameters A zu beachten. Die Arie „Vogelf£inger" 
10 ist mit einem mittleren Signalpegel von -32,15 dB auf Audio-CD gespeichert 
und bietet aufgrund von Signaldynamik und Klangfarben (Vorspiel, Gesang, 
Rohrfl5te) ein anspruchsvolles Beispiel zur Audiosignalcodierung. Da hier ein 
universelles PrddiktorTilter niedriger Ordnung (P « 2) verwendet wurde, 
kfinnen diese Simulationsergebnisse als reprfisentativ fOr eine Vielzahl von 
15 Audiosignalen betrachtet warden. 

Fig. 6 zeigt die gemessenen Stttrabstande, gemittelt Qber das ganze 
MusikstOck. Aus einem Vergleich der Fig. 3 (oder (1)) und Fig. 6 bei O = 1, /? 
= 7 ist offensichtlich, dass der mittlere Prddiktionsgewinn fOr diesen einfachen 
Pradiktor etwa 20 dB bis 23 dB betrdgt (Abtastfrequenz: 44,1 kHzl). Bei R'7 
20 wird der Gesamtgewinn bei alien Werten von D gut durch die Summe der 
beiden Einzelgewinne aus sphdrisch logarithmischer Quantisierung und 
DPCM approximiert, wobei dies bei R = 4 nur fUr D ^ 5 gilt. (Selbstverstdndlich 
ist (1) nicht anwendbar fUr beispielsweise /? - 3 Oder und D = 1, A = 
48270.) 

25 Der mittlere Signalpegel der OuvertOre liegt bei -27,20 dB; dieses MusikstOck 
ist durch einen sehr groflen Dynamikbereich von unterhalb -70 dB bis -17 dB 
charakterlsiert, siehe auch Fig. 7. In Fig. 7 ist der segmentweise Signalpegel 
sowie die SNRs fOr das Beispiel OuvertOre aufgetragen; jedes Segment 
besteht aus 6000 Abtastwerten und entspricht somit einem Zeitbereich von 

30 0,136 s. Hier wurde sphSrisch logarithmische Quantisierung in O = 3 
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Dimensionen mit A =1014 und A = 4858 sowie in O = 6 Dimensionen mit A 
s 48270 und dem gleichen Kompromiss-PrMiktor der Ordnung P » 2 wie Im 
ersten Beispiel angewendet. Die obere Kun^e zeigt die gro&e Dynamite, 
hervorgemfen durch die belcannten tauten T6ne getrennt durch 

5 Generalpausen Qber metirere Talcte zu Beginn und in der Mitte dieses 
popuiaren IS/lusitcstOclcs. Die Kurven fOr i4 = 1014 und A = 4858 offenbaren 
den Vorteil eines waclisenden A bzgl. des segmentweisen SNR, man 
beachte hier insbesondere die Generalpausen. Die untere Kurve in Fig. 7 
demonstriert diesen Aspekt fOr eine weitere Erli6hung von A und dem damit 

10 verbundenen weiteren SNR-Gewinn von 3 dB aufgrund der liOheren 
Dimensionszahl D = 6 (sielie auch Fig. 3). Es gilt zu beachten, dass trotz der 
niedrigen Rate von R - 4 bit/Abtastwert ein 10 logto(SNR) >35dB (bzgi. der 
originalen CD) selbst in den Generalpausen aufrecht erhatten wird. 

Der Wert A - 48270 entspricht einer Komprimlerung um 12 bits, z.B. fQr O = 
15 1, R 12 wird eine AuflSsung von 24 bits fQr kteine Signalwerte erzielt. 
Selbstverstandlich ist fOr D - 1, R = 4 wegen der zu geringen Intervatlanzahl 
eine Auflfisung entsprechend 16 bit niemals vorlianden. Aber bei Anwendung 
sphSriscli togarithmischer Quantislerung in 6 Dimensionen ist der extrem 
grofie Wert von A = 48270 fOr einen liohen Minlmalwert des SNR tatsdchlich 
20 gut gewahlt. Fur diesen Fall ist die Aufl5sung sogar identiscli mit derjenlgen 
der originalen CD-Daten fOr Segmente mit einem Signalpegel von -70 dB und 
niedriger. Bisher wurde selbst unter trainierten Personen niemand gefunden, 
der in der Lage war, zuverlSssig einen Unterscliied zwischen dem 
quantisierten Signal bei R = 4, D - 6, A = 48270 und dem Originalsignal zu 
25 hdren. 

Vorstehend wurde ein wellenfomierhaltendes Digitalisierungsverfahren fQr 
analoge Quellensignale gezeigt, das zum einen den Gewinn durch 
mehrdimensionale Quantisierung mit den Vorteilen der logaritlimischen 
Quantislerung verknQpft und zum anderen in der Lage ist, den objektiven 
30 Stdrabstand durcti Addition von Prddiktionsgewinnen zusatzlicii zu erhfitien. 
AbschlieKend sei angemerkt, dass sich dieses Verfahren neben einem 
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gOnstigen Austausch von Rate und Verzerrung vor allem durch einen extrem 
groften Dynamikbereich bei einer gleichzeitig Siu&erst geringen 
systembedingten VerzOgerung des Signals um nur wenige Abtastintervalle 
auszeichnet. 
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AnsorOche 

1 . Verfahren zur Digitalisierung analoger Quellensignale. mit dem Schritt: 

• DurchfOhren einer O-dimensionalen sphSirisch logarithmischen 
Quantlsierung der analogen Quellensignale. 

5 

2. Verfahren nach Anspruch 1 , wobei 

eine Transformation der Quellensignale in Kugelkoordinaten 
durchgefOhrt wird. 

10 3. Verfahren nach Anspruch 2, wobei 

eIne logarithmische Quantlsierung des Betrages in Kugelkoordinaten 
erfolgt. 

4. Verfahren nach Anspruch 2 Oder 3, wobei 

15 eine ungleichrndftige Quantlsierung der WInkelkoordinaten der 

Kugelkoordinaten durchgefUhrt wird. 

5. Verfahren nach Anspruch 1, mit den Schritten: 

- Quantisieren der analogen Quellensignale In Polarkoordinaten, 
20 wobei der Betrag der Quellensignale in Polarkoordinaten 

logarithmisch quantisiert wird. und wobei die Quantlsierung des 
Winkels In Polarkoordinaten der analogen Quellensignale In 
Abhangigkeit der quantislerten Winkel hOherer Ordnung erfolgt. 

25 6. Verfahren nach Anspruch 5, wobei 

die Quantisierung einer £>-dimensionalen Kugei mit einem 
Einheitsradius durch ein Netz aus D - 1-dimensionalen WQrfeln 
erfolgt. 



30 



Verfahren nach Anspruch S Oder 6. wobei 

diejenigen Quantisierungszellen ausgewShlt werden. welche auf der 
Oberfiache einer Kugel mit Einheitsradius liegen, wobei die 
Quantisierungszeilen O- 1-dimensionaie WUrfel darsteiien. 

Verfahren zur Komprimierung analoger Signale, mit den Schritten: 

- DurchfQhren einer Digitalisierung analoger Quellensignale gemdH 
einem Verfahren nach einem der AnsprQche 1 bis 7, und 

- DurchfQhren einer differenziellen Pulscodemodulation. 

Verfahren nach Anspruch 8, mit den Schritten 

- DurchfQhren einer Vorwdrtsprddii^tion zur Ermittiung eines 
Startwertes fQr Abtastwerte der Quantisierung auf der Basis des 
ai^tuellen Zustandes eines PrSdil^atorfilters, 

- Relconstruktion der nachfolgenden O Abtastwerte, 

- DurchfQhren einer 0-dimensionalen logarlthmischen sphSrischen 
Quantisierung der durch die VonA/Srtspradil^tion erhaltenen Werte, 
um eine Startzelle zu bestimmen, 

- iteratives Durchlaufen der Prddilction der differenzieilen 
Pulscodemodulation, um eine Quantisierungszelle mit dem 
geringsten Quantisierungsfehler zu bestimmen. 

Verfahren nach Anspruch 8 oder 9, wobei 

eine Lattice-Decodierung zur Bestimmung einer gQnstigen 
Quantisierungszelle durchgefQhrt wird. 

Verfahren nach einem der AnsprQche 8 bis 10, wobei 

eine RQckwdrtspradilction einer differentiellen Pulscodemodulation 
durchgefQhrt wird, wobei nach der Verarbeitung von D Abtastwerten 
eine Aktualisierung des Pradiktionsfiiters um D Schritte erfolgt. 
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12. Vorrichtung zur Digitalisieaing analoger Quellensignale, mit 

- einer 0-dimensionalen sphSrisch logarithmischen 
Quantisienjngseinrichtung. 

5 

13. Vorrichtung nach Anspruch 12, mit 

einer Transformationseinrichtung zur DurciifUhrung einer 
Transformation der analogen Quellensignaie in Kugetl^oordinaten. 



10 14. Vorrichtung nach Anspruch 12 Oder 13, wobei 

eine logarithmische Quantisierung des Betrages erfolgt. 

15. Vorrichtung nach Anspruch 13 Oder 14, wobei 

eine ungleichmSl^ige Quantisierung der Winlcellcoordinaten der 
15 sphSrischen Oder Kugellcoordinaten durchgefOhrt wird. 



16. Vonichtung zur Komprimierung analoger Signale mit einer Vorrichtung 

zur Digitalisierung analoger Signale nach einem der AnsprQche 12 bis 
15, und einem Encoder fOr differenzielle Pulscodemodulation. 

20 

17. Vonichtung nach Anspruch 16, mit 

- einer Vonvartspradil<tionseinrichtung zur Ermittlung eines Startwertes 
fUr die Abtastwerte der Quantisierung auf der Basis des alctuellen 
Zustandes eines Pradil<atorfilters, 

25 - einer Rekonstrulctionseinrichtung zur Rekonstruktion der 

nachfolgenden D Abtastwerte, 

- einer O-dimensionalen logarithmischen sphSrischen 
Quantisierungseinrichtung zur Quantisierung der durch die 
VorwSrtsprddiktion erhaltenen Werte, um eine Startzelle zu 

30 bestimmen, 

wobei die PrSdiktion der differenzielien Pulscodemodulation iterativ 
durchlaufen wird, um eine Quantisierungszelle mit dem geringsten 
Quantisierungsfehler zu bestimmen. 
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1 8. Vorrichtung nach Anspruch 1 7 Oder 1 8, wobel 

eine Lattice-Decodierungseinrichtung zur Bestimmung einer 
gOnstigen Quantisierungszelle durchgefOhrt wird. 

1 9. Vorrichtung nach einenrt der AnsprUche 1 7 bis 1 9, wobei 

eine RQckw^rtsprddiktion einer differentielten Puiscodemodulation 
durchgefOhrt wird, wobei nach der Verarbeitung von D Abtastwerten 
eine Aktualisierung des Prfldil(tionsfiiters um O Schritte erfolgt. 
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Fig. 2 
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Fig. 3 
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Fig, 5 
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